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ABSTRACT 

We determine photometric metal abundance estimates for individual main-sequence stars in the Virgo 
Overdensity (VOD), which covers almost 1000 deg 2 on the sky, based on a calibration of the metallicity 
sensitivity of stellar isochrones in the gri filter passbands using field stars with well-determined spectro- 
scopic metal abundances. Despite the low precision of the method for individual stars, we derive [Fe/H] = 
-2.0 ±0.1 (internal) ±0.5 (systematic) for the metal abundance of the VOD from photometric measurements 
of 0.7 million stars in the Northern Galactic hemisphere with heliocentric distances from ~ 10 kpc to ~ 20 kpc. 
The metallicity of the VOD is indistinguishable, within A[Fe/H] < 0.2, from that of field halo stars covering 
the same distance range. This initial application suggests that the SDSS gri passbands can be used to probe 
the properties of main-sequence stars beyond ~ 10 kpc, complementing studies of nearby stars from more 
metallicity-sensitive color indices that involve the u passband. 

Subject headings: Galaxy: abundances — Galaxy: evolution — Galaxy: formation — Galaxy: halo — Galaxy: 
stellar content — Galaxy: structure 



1. INTRODUCTION 

The structure, chemistry, and kinematics of the stellar halo 
of the Galaxy, with its predominantly old and metal-poor 
populations, collectively preserve a detailed re cord of our 
Galaxy's formation in the early universe (e.g., E ggen et"aLl 
[19621 LSearle & Zinnll978l) . Thanks to larg e -area surveys such 
as the Sloan D i gital S ky Survey (SPSS: lYorketall 1 20001: 
Stoughton etail 120021: lA bazaiian etap 120031 12004 l200l 
Adel man-McCarthv et al.ll2006l 120071 l2008t lAbazaiian et all 
2009), recent studies have revealed that the halo is marked by 
numerous stellar substructures. The presence of these lumpy 
and complex substructures (both in real and phase space) are 
in qualitative agreement with models for the formation of the 
stellar halo through the hierarchical merging and ac cretion of 
low-mass subhalos (e.g., Bul lock & Joh nston 20051). 

Among the various substructures discovered to date, the 
Virgo Overdensity (VOD) is one of the most striking. It was 
discovered as a stellar overdensity of main-sequence stars in 
SDSS; starcounts in the region toward Virgo are enhanced 
by a factor of two above the background stellar distribution 
dJuric et al.ll2008l) . The overdensity covers almost 1000 deg 2 , 
and it appears to span a wide range of heliocentric distances 
of ~ 5 kpc-20 kpc. The over density seems to be associated 
with c lumps of RR Lyrae stars (Vivas et al. 2001; Duff au et al.l 
120061) and turn-off stars dNewberg et al.l 120021) . but it is less 
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likely to be connecte d with the leading t idal tail of the Sagit- 
tarius dwarf galaxy (Newberg et al. 20Q7[). 

At present, metallicity estimation from broadband photom- 
etry is the only practical means of obtaining metal abundances 
for a large number of faint objects such as those in the VOD. 
Such methods are based on the relative sensitivity of stellar 
colors to photospheric abundances over a wide wavelength 
baseline. The clear advantage of using a photometric metal- 
licity technique is the efficiency of estimating metallicities for 
individual main-sequence stars, which are the most plentiful 
a nd representative sam ple of stellar populations. 

Ilvezic et ail d2008al) constructed photometric metallicity re- 
lations in the u-g vs. g—r plane using SDSS filter passbands, 
and studied the abundance structures of the Galaxy with an ac- 
curacy of ~ 0.2 dex at g < 17. This approac h is similar to the 
traditional UBV method (e.g., Carnev 1979), which relies on 
the strong dependence of f/-band magnitudes on metal abun- 
dance. However, the M-band photometry in SDSS is limited 
to u rj 22 (99% detection limit). This, and the greatly deterio- 
rating errors in SDSS w-band magnitudes near the faint limit, 
restricts photometric metallicity estimates to stars at r < 20.8, 
an insufficient d epth to fully explore the VOD (see Fig. 37 in 
Uuric et al.ll2008l) . 

In this letter we overcome the limitations of the M-band pho- 
tometry in SDSS by exploring less metallicity-sensitive, but 
better-determined, color indices in the gri passbands. 9 Turn- 
off stars in globular clusters have M r ~ 4, so the SDSS survey 
limit in gri (95% completeness limit at r = 22.1) allows us to 
probe the halo out to ~ 20 kpc using stars that are ~ 2 mag 
below the main-sequence turn-off (~ 0.6 M Q ). 

2. METHOD 

2.1. Photometry 

We employed SDSS photometry from DR7 
(AbazajianetalJ[200l). SDSS measures the brightnesses of 

9 The z-band data were not used due to the bright survey limit in SDSS for 
this filter. 
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stars with accurate astrometric po sitions dPier et alj [2003). 
using a dedicated 2.5- m telescope (|Gunn et alj 12006ft in five 
broadband filters (ugriz; Stoughton et al. 2002), on 6 columns 
of CCDs (iGunn et al] 119981) . under photometric conditions 
(Hogg et alj 120011) . Photometric calibration is carried out 
using observations of sta rs in the secon dary patch transfer 
fields dTucker et al]|2006l) , based on the ISmith etail (120021) 
sample of standard stars. 

The rms photometric precision is 0.02 mag for sources not 
limited by photon statistics, and the photometric calibration 
is accurate to ~ 2% in the gri bands, a nd ~ 3% in u and z 
dlvezic et alj|200l l2004t lAn et alj|2008h . We used photome- 
try of stellar objects (i dentified as S TAR in the standard SDSS 
photometric pipelines; Lupton et al. 2002) that were detected 
(at the 5cr level) in all of the gri passbands. The observed 
magnitudes are corrected for exti nction adopting reddening 
values in the Schleg el et alJ (11998ft dust maps and the extinc- 
tion coefficients given by lAn et alJ d2009h . 

2.2. Photometric Metallicity 

The isochrones in lAn et alj d2009ft were used to deter- 
mine photometric metallicity ([Fe/H] p h t) estimates based 
on color-color relations for main-sequence stars. In the 
following analysis, we a dopted th e same a-element en- 
hancement scheme as in lAn et alj (|2009), motivated by 
the observed behavior of these elements among field 
dwarfs and clust e r stars from spectro scopic studies (e.g., 
iVenn et alJ l200fl iKirbv et all l2008ah : [o/Fe]= +0.4 at 
[Fe/H]= -3.0, [a/Fe]= +0.3 at [Fe/H]= -2.0,-1.5,-1.0, 
[a/Fe]= +0.2 at [Fe/H]= -0.5, and [a/Fe]= +0.0 at [Fe/H]= 
-0.3,-0.2,-0.1,+0.0,+0.1,+0.2,+0.4. A linear interpolation 
was used in this metallicity grid to obtain isochrones at inter- 
mediate [Fe/H] values. We adopted an age of 12.6 Gyr for 
[Fe/H] < -1, and 4.6 Gyr at [Fe/H] > 0, and a linear inter- 
polation between these two values. To expedite the computa- 
tional process, we derived 5 ,/! -order polynomials to describe 
the color-magnitude relations, with intervals in abundance of 
A[Fe/H] = 0.01 dex. 

To reliably estimate a photometric metal abundance, it 
is necessary to use stellar models that match cluster main- 
sequences and give internally consistent distances from mul- 
tiple color-magnitude diagrams (CMDs). However, despite 
improvements in theoretical models, there still exist small 
but significant mismatches between calculated isochrone col- 
ors and the best photometry in well-stud ied open clusters 
dPinsonneault et al.ll2004t lAn et al.ll2007allbb . Therefore, we 
applied empirical corrections on theoretical color-!T e ff rela- 
tions for [Fe/H] > -0.8 to match the observed main se quence 
of the solar-metallicity cluster M67 dAn et al.l 12009ft . Con- 
stant correction factors were adopted at [Fe/H] > 0, and a 
linear ramp was used between [Fe/H] = -0.8 and [Fe/H] = 0, 
so that the correction becomes zero at [Fe/H] = -0.8 and be- 
low, as the mode ls are in good agreement with the data for 
globular clusters dAn et 
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A photometric metallicity (lAn et a l. 2007b) was computed 
for each star by requiring distances from main-sequence fit- 
ting to be the same from two different CMDs, with g — r and 
g — i as color indices and r as a luminosity index (Fig. [TJ. 
We searched the entire [Fe/H] grid from -3.0 to +0.4 to find 
[Fe/H] with a minimum \ 2 value, defined as 
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FIG. 1. — Illustration of the photometric metallicity estimation technique. 
Top: Distance moduli from two CMDs as a function of a metallicity. For 
this example, [Fe/H] = —1.6 was assumed for the true metallicity of a star. 
Bottom: Difference in distance modulus from the two CMDs. A photometric 
metallicity is defined as the [Fe/H] that results in the same distance modulus 
from the two CMDs. The gray region shows the error bound for a ~ 2% 
photometric color error. Note that only the upper [Fe/H] bound is defined in 
this example. 

for each star. Here, /i and tr M are the distance modulus and 
its error for each CMD, respectively. The quantity j2 is a 
weighted average distance modulus from the (g - r, r) and 
(g-i,r) CMDs. Since only three passbands are considered 
here, the problem is reduced to the traditional manner of de- 
termining metal abundances from a color-color diagram. In 
some cases, a star becomes bluer than the main-sequence 
turnoff as the metallicity increases, either due to a large pho- 
tometric error or a younger age than our assumed values in 
the models. If a minimum % 2 was not found, [Fe/H] p h ot 
was estimated based on a mean relation between [Fe/H] and 
A(m-M) . 

2.3. Accuracy of Photometric Metallicity Determinations 

Photometric metallicity estimates become insensitive for 
very metal-poor stars, as illustrated in Figure Q] Therefore, 
[Fe/H] p hot could be biased due to systematic effects such as 
photometric errors, unresolved binary stars, and dust extinc- 
tion, which alter the differential color estimates at a small 
level. We tested our photometric metallicity estimates using 
field stars with well-measured spectroscopic estimates, as dis- 
cussed below, and adjusted [Fe/H] p h t to correct for the bias. 

Figure [2] shows the comparison with low-resolution 
(R ~ 2000) spectroscopic measurem ents for field dwarfs 
dCarollo et al.ll2007t lYannv et alll2009ft from the most recent 
version (Data Release planned in December 2010) o f SSPP 
(SEGUE Stellar Parameter Pipeline; lLee et alj|2008albb . This 
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FIG. 2. — Comparison between spectroscopic and photometric estimates 
of [Fe/H]. Open circles are median [Fe/H] p hot for the SSPP field star sam- 
ple with available spectroscopic metallicities; the gray region represents in- 
terquartile ranges. Solid lines are piecewise linear fits to the open circles. 



version of the SSPP partially solves the under- and over- 
estimation of metallicity at higher and lower ends of [Fe/H] 
that the earlier version (DR7) showed. Comparisons are 
shown for 46,983 stars, after applying cuts at S/N > 20 and 
logg > 4.2 to select dwarfs with good spectroscopic abun- 
dance measurements. We further restricted the sample to 
those with |1.4(g — r)o — (g — Oo| < 0.1. Open points show a 
median metallicity in [Fe/H] = 0.2 dex bins; the gray region 
represents interquartile ranges. Because the metallicity sen- 
sitivity essentially disappears below the lower limit of our 
metallicity grid, a significant number of stars are found at 
[Fe/H] phot = -3.0. 

The large uncertainty in [Fe/H] p h t is mainly due to photo- 
metric errors. We performed artificial star tests by generating 
stars from the isochrones with Gaussian errors of 0.02 mag in 
r, g-r, and g—i. The interquartile range from the test showed 
a reasonable agreement with those for the SSPP sample, in- 
dicating that the observed dispersion is at least in part due to 
the ~ 2% photometric error in SDSS. Despite the large un- 
certainty in [Fe/H] p hot for individual stars, a meaningful con- 
straint on a median [Fe/H] p h t can be made by applying the 
technique to a large number of stars (N ~ 500, see below). 

For the SSPP sample in Figure|2] we estimated the error in 
the median of < 0.1 dex at [Fe/H] spec > -2 in each [Fe/H] 
bin; this was done by computing the median absolute de- 
viation (MAD) 10 from stars with [Fe/H] p h t above the me- 
dian, divided by y/N (although [Fe/H] p h t does not strictly 
follow a normal distribution). We explored the effects of dif- 
fering stellar ages (a = 20%), [a/Fe] ratios (a w 0.1 dex at 
[Fe/H] < — 1), and dust extinctions (a = 20%), but they change 
[Fe/H] phot by<0.1dex. 

As shown in Figure [2] our photometric solution underesti- 
mates [Fe/H] by as much as A [Fe/H] ~ 0.5 dex for [Fe/H] < 

10 MAD = 1.483 median(|.v, -median(x,)|). 



-1 with respect to the SSPP results. One reason for this off- 
set is likely due to the presence of unresolved binary pop- 
ulations in the sample, which will have a p ortion of their 
light from a cooler secondary (lAn et alJl2007bh . To evaluate 
the effect of binaries, we performed artificial star tests with 
binaries generated from a flat mass function for the secon- 
daries, and found that the A[Fe/H] ~ 0.5 dex offset can be 
explained with a ~ 40% binary fraction. On the other hand, 
the difference could also be induced by a systematic color 
offset in the isochrones. Although we found a good agree- 
ment ofthemodels with globular cluster data within the er- 
rors dAn et alj 12009). a small color offset (< 0.01 mag) in 
the models can still change the photometric metallicity esti- 
mate by A [Fe/H] ~ 0.5 dex at the low metallicity end. Al- 
though possible, the low-resolution spectroscopic values are 
less likely the source of the problem, given the extensive 
tests th at have been applied in their validation (see iLee et alj 
I2008allbl ; lAUende Prieto et al.| [2008). Nevertheless, the rela- 
tive metallicity comparison would be robust for stellar popu- 
lations with the same binary fraction, even if there is an offset 
with the SSPP. 

In the following analysis, we adjusted our [Fe/H] p h t to be 
on the same scale as the SSPP results by deriving a piecewise 
linear fit to the median [Fe/H] p h t as a function of metallic- 
ity (solid lines in Fig. [2]). This adjustment implicitly assumes 
the same fraction of unresolved binaries applies to both the 
SSPP and the field halo samples, if binaries are solely respon- 
sible for the difference between the [Fe/H] p h t and the SSPP 
results. No statistically significant trend was found in the sys- 
tematic difference between [Fe/H] p h t and [Fe/H] spec for var- 
ious sets of the binned data with different magnitude and color 
ranges. 

3. RESULTS 

Figure [3] shows the distribution of stellar number density 
(top) and photometric metallicity (bottom) for 740,658 stars 
detected in SDSS. The maps are Lambert equal-area projec- 
tions of the Northern Galactic hemisphere, and the distribu- 
tions are projected as seen from the Sun (i.e., a view from 
the inside of the Galaxy). The North Galactic Pole is located 
at the center, and the direction of the Galactic Center is to- 
ward the bottom in each panel. Each pixel has an area of 
12.96 deg 2 . To avoid any possible bias, we restricted our sam- 
ple to stars with 0.3 < (g— r)o < 0.4 and the same color-color 
cut as those for the SSPP comparison sample. We used stars 
with 15.0 < (m-M)o < 16.5 (10 <d< 20 kpc). Although a 
distance modulus was derived from the \ 2 minimization for 
each star, we used a 12.6 Gyr old model with [Fe/H] = -1.6 
for all of the stars to bracket the distance range in the sample. 
This was because the uncertainty in [Fe/H] p h t for individ- 
ual stars is large enough that the strong correlation between 
[Fe/H]ph ot and distance could lead to a biased resu lt. 

The feature at (l,b) ~ (300°, 70°) is the VOD dJuric et alj 
2008), and the long strip that crosses the sky from (l,b) ~ 
(210°, 30°) to the VOD is the leading tidal tail of Sgr. The 
feature at (l,b) ~ (40°, 40°) is the Hercules-Aquila Cloud, 
anothe r large-area overdensit y of halo stars discovered from 
SDSS (Bel okurov et alj|2007l) . 

The median metallicity of stars in the VOD area at 270° < 
/ < 330° and 60° < b < 70° is [Fe/H] = -2.0 ±0.1 from the 
metallicity map in Figure[3] where the error is from a pixel-to- 
pixel dispersion. The photometric zero points in SP SS vary at 
the ~ 2% level over a ~ 10' scale along the scan line dAn et alj 
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FIG. 3. — Top: Number density of stars at distances from ~ 10 kpc to 
~ 20 kpc from the Sun in the Lambert projection of the Galactic coordinates. 
The North Galactic pole is at the center, and the Galactic Center is to the 
bottom. Concentric circles represent b = 0° , 30° , and 60° , respectively. The 
VOD is the feature seen at (l,b) ~ (300°, 70°). Bottom: Median metallicity 
of the same stars as in the top panel. The median occupancy of each pixel is 
544 for both maps. 

2008). They are also known to vary at the same ~ 2% level 
over a larger angular scale along the stripe (a 2.5° wide SDSS 
stripe typically runs from the I s t to the 3 rd Galactic quad rant), 
and from one strip to the other dlvezic et alj |2003, 2004). Our 
hope is that these components are properly averaged out over 
a large area in the sky, such as the solid angle covered by the 
VOD. 



TABLE 1 

Error budget in [Fe/H] phot for the VOD 



Source of Error 


A Quantity 


A[Fe/H] phol 


Internal 




±0.1 


Age 


±20% 


±0.3 


Dust extinction 


±20% 


±0.1 


[ce/Fe] 


±0.1 dex 


±0.3 


Contamination by giants 




-0.1 


Alternative approach 




-0.1 


Total (systematic) 




±0.5 



NOTE. — These estimates are slightly different from those for the SSPP 
sample (§ 2.3) due to different color ranges used. 



Field halo stars located at the mirrored position (30° < / < 
90°, 60° < b < 70°) exhibit [Fe/H] = -1.9 ±0.1. Although 
half of the stars in the direction toward Virgo are likely as- 
sociated with a progenitor dwarf galaxy or a tidal stream, the 
median [Fe/H] value essentially remains unchanged; it is the 
same as that for the field halo stars within the precision of the 
technique. 

It is perhaps of interest that the field-star metall icities are 
as low as they appear to be, as ICarollo et alJ (120071) have ar- 
gued that the peak metallicity of the outer-halo population is 
[Fe/H] = -2.2, and that this component is expected to dom- 
inate over the more metal-rich inner-halo population (with 
a peak metallicity at [Fe/H] = —1.6) at Galactocentric dis- 
tances greater than 15 kpc - 20 kpc. Note that metallicity 
estimates at the second Galactic quadrant are even lower than 
[Fe/H] « -2. 

In Table Q] we list the sources of several systematic errors 
and their contributions to errors in [Fe/H] p h t for the VOD. 
The effects of the age, extinction, and [a/Fe] were tested 
by constructing similar [Fe/H] p h t maps to Figure[3]with dif- 
ferent parameters in the models. Our photometric technique 
is valid only for main-sequence stars, but giants constitute 
appro ximately 10% of the stars in the sample (see lAn et alJ 
2008). Since photometry alone cannot be used to adequately 
distinguish giants from dwarfs, we estimated a bias due to the 
presence of giants using photome try for a sample of nearby 
globular clusters dAn et al.1 120081) . In Table Q] we also list 
the error from an alternative approach where we used a me- 
dian difference in distance modulus in each pixel of Figure [3] 
to determine [Fe/H] p h t- The total systematic error is the 
quadrature sum of all of the error contributions. The rel- 
ative metallicity comparison is more robust, if stellar pop- 
ulations in the field halo and the VOD have the same age 
and [a/Fe] distributions: the difference in [Fe/H] p h t remains 
within A[Fe/H] phot < 0.2. 

4. DISCUSSION 

This initial application of the photometric metallicity tech- 
nique demonstrated that there is no metallicity difference be- 
tween the field halo stars and those in the VOD within the 
precision of the method. Our estimate can be compared 
with previous measurements for a handful of RR Lyrae vari- 
ables that are likely associated with the V OD dDuffau et alJ 
2006: 1 Vivas etani200l IPrior et alJl2009l ). These estimates 
range from [Fe/H] = -1.55 to -1.95, based on the pseudo- 
equivalent width of the Ca II K line. From p rincipal axes on 
the u—g vs. g—r diagram. lJuric et aTl {2008) argued that the 
VOD metallicity is lower than that of thick-disk stars and sim- 
ilar to that of halo stars. 

It is tempting to place the VOD in the observed trend of the 
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luminosity -metallicity relati on among surviving dwarf galax- 
ies in the Local Group (e .g., iGrebel et alj|2003t iMufioz et al.l 
I2006t iKirbv et alj2008bb. If we take our median [Fe /HI from 
Figure E] at M v ~ 10 (lJuric et"aTll2008t IPrior et al.ll2009l) . the 
VOD follows this trend, supporting the idea that [Fe/H] can 
serve as a luminosity indicator for an accreting dwarf galaxy, 
in the process of building up the stell ar halos of large spi - 
ral galaxies like the Milky Way (e.g.. iJohnston et ail 2008). 
However, this should be taken with caution, as previous stud- 
ies often report average abundances for dwarf galaxies rather 
than median values. 

Future imaging survey s, such as the the L arge Synoptic 
Survey Telescope (LSST; ITvezic et al]|2008bl) will use simi- 
lar photometric bandpasses as those in SDSS, providing even 
deeper (and far more accurate) photometric data than SDSS 
over a larger fraction of the sky. Our photometric metallicity 
method will be useful to exploit these databases for under- 
standing the chemical evolution of progenitor dwarf galaxies 
that are identified, as well as for the bulk populations of field 
stars. 
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